Holographic Operator Mixing and 
Quasinormal Modes on the Brane 



Matthias Kaminski a \ Karl Landsteiner 6 , Javier Mas c , Jonathan P. Shock c and 
Javier Tarn'o c 

a Department of Physics, Princeton University 
Princeton, NJ 08544, USA 

b Instituto de Fisica Teorica CSIC-UAM 
C-XVI Universidad Autonoma de Madrid 
E- 28049 Madrid, Spain 

c Depto. de Fisica de Particulas, Universidade de Santiago de Compostela, and 
Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E- 15782 Santiago de Compostela, Spain 

e-mail: mkaminsk@princeton.edu , karl.landsteiner@uam.es, 
jamas@fpaxpl.usc.es, shock@fpaxpl.usc.es, tarrio@fpaxpl.usc.es 
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from general bilinear actions and fields obeying coupled differential equations in the 
bulk. The matrix-valued spectral function is shown to be independent of the radial 
bulk coordinate. Applying this framework we improve the analysis of fluctuations 
in the D3/D7 system at finite baryon density, where the longitudinal perturbations 
of the world-volume gauge field couple to the scalar fluctuations of the brane em- 
bedding. We compute the spectral function and show how its properties are related 
to the quasinormal mode spectrum. We study the crossover from the hydrodynamic 
diffusive to the reactive regime and the movement of quasinormal modes as functions 
of temperature and density. We also compute their dispersion relations and find that 
they asymptote to the lightcone for large momenta. 
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1. Introduction 



The regime of application of the gauge/gravity correspondence [1] is quickly being 
extended to cover ever more diverse areas of theoretical physics. Its utility derives 
from the fact that it is a weak/strong coupling duality, meaning that by studying 
classical gravity we can obtain information about strongly coupled gauge theories. 
One reason for the growth in this applicability is due to the fact that the correspon- 
dence allows the calculation of real time Green's functions at finite temperature and 
density [2]. This paves the way to obtaining physical transport coefficients, such as 
the shear viscosity and the conductivity. 

The original conjecture relating M = 4 Super Yang Mills to type IIB super- 
gravity in AdS^ x S 5 has been generalized to cover many other gauge theories by 
considering different gravitational backgrounds (not necessarily derived from critical 
string theory, e.g. [3]). Usually, such generalizations only include fields in the adjoint 
representation of the gauge group. It was shown in [4] how to add matter in the fun- 
damental representation, by considering D7-branes extended along the Minkowski 
and radial directions and wrapping an S 3 inside the S 5 . These fundamental branes 
are typically studied in the limit in which they do not backreact on the geometry 
of the D3-branes. This means that the quarks do not contribute to the dynam- 
ics beyond tree-level - it is the quenched approximation. Further generalizations 
added temperature to the system by considering black _D3-branes, holographically 
describing a thermal field theory [5]. 

The finite temperature plasma phase of strongly coupled M = 4 SYM is modelled 
as an asymptotically AdS black hole with planar horizon topology. The D7-branes 
can be embedded in two rather different ways with completely different gauge theory 
phenomenology. If the D7-brane does not enter the horizon the fields on it have 
a discrete spectrum of normal modes whose dual interpretation is as stable mesons 
(quark-antiquark bound states - nb. not confined) [6]. Instead , if the D7-brane 
reaches the horizon, the modes are quasinormal with complex frequncies [7]. In 
certain ranges of the parameter space [8] these quasinormal modes lie close to the 
real axis. Then the spectral function exhibits sharp peaks that we may interpret as 
quasiparticle states. These 'black hole' embeddings are then naturally associated to 
the high temperature phase (with respect to the quark mass) in which the quarks 
can not form long lived bound states - the mesons melt into the surrounding adjoint 
plasma. 

A particularly interesting generalization of this setup is provided by the inclusion 
of finite baryon density [9, 10]. Due to the presence of electric flux on the D7-brane 
- the holographic dual of finite baryon density - the brane is forced to end on the 
horizon. The flux can not end on a shrinking S 3 and has to fall into the black 
hole. The study of the meson spectral functions in this setting has attracted interest 
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in the past and it was shown that for low densities and low temperature to quark 
mass ratios the spectral functions present peaks corresponding to long lived states 
that approach the zero temperature meson masses of the Minkowski embeddings at 
vanishing density [11]. 

As has been appreciated before [12, 13], the fluctuations of the longitudinal 
vector and scalar components couple when the effects of baryon density are taken 
into account, making the analysis of this sector much more involved. Physically the 
coupling occurs because a perturbation in the scalar sector is a modification of the 
brane embedding and this has to backreact on the charge distribution on the brane, 
which in turn modifies the form of the fluctuation equations for the electric fields on 
the brane. 

At finite density all consistent embeddings are of the black hole type and the 
dynamics of the system can be described in terms of its quasinormal modes. Since 
we now have to deal with a system of coupled differential equations, we ought to 
state precisely the new prescription that generalizes the computation of quasinormal 
modes. A similar problem has been considered before in the context of holographic 
superconductors, and here we extend the formalism put forward in [14] to a generic 
bilinear bulk action. We explictely construct the matrix of retarded Green's func- 
tions. From it, the spectral function matrix can be derived and we prove that it is 
fully independent of the radial bulk coordinate, a result that was known for single 
field correlators. In turn, the quasinormal modes can then be defined as zeroes of 
the determinant of a maximal set of linearly independent solutions evaluated at the 
boundary of AdS space. 

The paper is organized as follows. In section 2 we develop the formalism which 
allows us to define a matrix of retarded Green's functions for a coupled system from 
the original bilinear action. We show that the matrix-valued spectral function can 
be interpreted as the matrix of Noether currents of a collection of global U(l) bulk 
symmetries. Finally we explain how to calculate the quantities of interest in the case 
in which the solutions have to be found using numerical methods. 

In section 3 we summarize the description of the quenched D3/D7 system in the 
presence of baryon density, focusing on the fluctuations of the longitudinal electric 
field and the embedding profile. 

In section 4 we apply our general formalism to this system and describe the 
results obtained, focusing on the effects of the mixing in the system. 

2. General formalism and methods 
2.1 Holographic operator mixing 

In the holographic context, operator mixing under the RG flow manifests itself as a 
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coupled system of differential equations. In some contexts the system of equations 
may be separable, but generically, when the mixing matrices between the derivative 
and non-derivative terms differ from each other, this is not so. In the present case 
we focus on non-separable systems where an analysis of the coupled set of differen- 
tial equations is necessary. Consequently, holographically obtaining the sources and 
expectation values of the mixed operators corresponds to finding a particular set of 
solutions to the coupled system of equations. With this aim in mind we consider 
here a general bilinear bulk action for N fields $ 7 , I E {1, • • • , N} 

S = J d d x J dz [d m ^ I A IJ {x,z)d n ^ J 1 mn + &B™(x,z)d m ® J + <£> I C IJ (x,z)$ J ] , 

(2.1) 

where m, n span the Minkowski and radial coordinates (x ~ x M , z). i.e., we consider 
that the fields &(x,z) and the matrices Au(x,z), Bjj{x,z) and Cjj(x,z) have no 
dependence on any transverse coordinate, which we integrate out in the action S. 
Furthermore, apart from being real, no symmetry properties will be assumed for 
these couplings 1 . Inserting the Fourier transform 



*V,z)= / j—.H (z)e-^, (2.2) 



d d k , 
(2vr) c 

into (2.1), standard manipulations lead to an action for the Fourier modes of the 
following general form 2 



S = J-^jdz [^\A u (k,z)^ + ^ k B u (k,z)^i + ^ k Cu(k,z)^ k 

(2.4) 

with k = Aij(— k,z) = Ajj(k,z)*, and equivalently for B and C. Now in 
order to avoid double counting, we split the momentum integration into "positive" 
(k > = (u;>0,q)) and "negative" (fc< = (u;<0, q)) momenta. Thus 

S = Jdhjdz [2Afj^' J k + BA$f + B\j#L k *l + 2Cfj$L k Q{\ , (2.5) 



^^This is because we are interested in a system with finite baryon density, modelled by the 
presence of a background Aq component of a U (1) gauge field. In such setups the 7 in equation 
(2.1). while not being symmetric, plays the role of the induced metric. 

2 In going from (2.1) to (2.4) we relate 

Au(k,z) = Af J ( a; ,z)7 zz , (2.3a) 
Bu{k, z) = ~2ik^ z Afj(x, z) + Bfj(x, z) , (2.3b) 
Cu(k, z) = -k^k v Ajj(x, z) - ik^Bfjix, z) + Cu(x, z) , (2.3c) 

where the superscript S(A) denotes (anti-)symmetrization i.e. Mfj A = | (Mjj ± Mji). Notice 
that in (2.5) we have admitted a slight generalization in which A has a k dependence. This will be 
the case when one performs complex valued changes of variables like the one for the gauge invariant 
combination En = ifwA 1 + qA°) which will be needed later. 
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where J dk> = J °° du J Rti _i <i d-1 q. Hereafter will always be assumed to be 
"positive", k = k > . Also M H,A now stands for the (anti-)hermitian part M H,A = 
\{M ± M^). Written in this form, a given mode, say $^ = n x 00 ), only enters once in 
each bilinear term 3 . Varying holding fixed, the Euler-Lagrange equations 
of motion follow 



[E.O.M. 



-2{Afj&i)' + 2Bfj& k J + (2C H - B%M = 



(2.6) 



Upon solving the equations of motion (2.6), one may find that asymptotically 



near the boundary 4 , the components of the vector $ go like $ 7 (z — > 0) 



z A ic 



!>{ + .... Namely A{_ is the smallest exponent at the boundary z = 0. In 
order to compute the Green's functions of the dual quantum operators we choose to 
consider conveniently normalized fields &k{z) — zA ~&l( z ) that close to the boundary 
have an expansion $ / (z — >• 0) = 0q + O (z a +~ a -) , meaning that </>q can be interpreted 
as the source of the dual operator 5 . The new fields can be treated collectively in the 
same formalism by defining the rescaling matrix D 1 j = 5 1 jz A - = D^ 1 j. Replacing 
$ by D<& inside (2.5) yields a new action of the same form 



S 



now with 



dk^ / dz 



Bjj¥_M j 



-k^k 



2C? J <5> I _ l .<$> J 



-k^k 



(2.7) 



A H 
B 

C H 



D ] A H D , 

D ] BD + 2D' ] A H D , 
D ] C H D + D' ] A H D' 



-D ] BD' + -D' ] B ] D . 



(2.8a) 
(2.8b) 

(2.8c) 



Hereafter we will assume without loss of generality that the fields are normalized in 
this way. For these normalized fields the action is given by (2.7) and henceforth we 
shall omit all bars. 



2.2 Holographic Green's functions 

We now want to construct the precise solutions which are sources for operators 
O 1 . The fact that the fields are solutions to a coupled system of differential equations 
can be interpreted as the holographic dual of operator mixing. This means in turn 
that we cannot simply speak of a single operator O 1 , but we must specify at which 

Alternatively one could use the reality condition $^ fe = and treat this field as independent 
of<&£. 

4 Notice that we choose the radial variable to present the boundary at z — > 0, thus the IR of the 
theory will be at a positive scale z%. 

5 Care must be taken to ensure that such a mode is in fact non-normalizablc. This is independent 
of the redefinition discussed above. 
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scale this is defined. The most natural choice is to define the operators in the UV 
at a cutoff scale za which will ultimately be taken to the boundary. In the example 
discussed here there will be a generic U(l) gauge symmetry present on the world- 
volume. Ultimately we want to think of the dual global U(l) symmetry as being 
weakly gauged such that the spectral functions of the conserved current can be used 
to calculate dilepton and photon production rates in a charged plasma along the lines 
of [15]. The correct operator to be coupled to the electromagnetic photon in the field 
theory is of course the one that sources only the electric current in the ultraviolet 6 . 

We now construct $^'s that are solutions to the coupled set of N differential 
equations in the bulk and whose boundary values serve as the sources for operators 
I (k). Concretely, let us set / = 1. A particular solution which sources O 1 will be 
given by a vector of functions (<& k (z), & k (z), ...) that as we approach the UV cutoff 
asymptotes to a single component vector, say (& k (z),$l(z), ...) (<^, 0, 0, ...). 

The same is true for any J = 2, 3, • • • , N. Hence, collectively, a bulk solution dual 
to a source O lo (k) is given by a set of functions {^j,(z)} which solves the equations 
of motion in the bulk and asymptote to ^ J (z\) = 5 J / o 0q°(/c), J = 1,...,N, where 
0q° (k) = ip 1 ^ is the source of the corresponding operator O lo (k). 

Because the system of differential equations is coupled, at any other scale z > 
this set of functions, {$£(z)}, will in general source a linear combination of all the 
operators. Hence, this set of functions can be written in terms of the boundary 
values, <p k , as follows 

& k (z) =F I J (k,z)<p J k , (2.9a) 
QUz) = F'j(-k, z)ifii k = ^/(k, z) , (2.9b) 

with (p k arbitrary (sourcing the corresponding operators) and all the dynamics of the 
fields encoded in the "solution matrix" F(k, z) 1 j = F(—k, z)* 1 j, normalized at the 
UV cutoff radius z&, such that 

F(k, z A y j — Sj. (2.10) 

Any complete set of independent solutions to the equations of motion is enough to 
build the matrix F, and we shall give a concrete prescription below. For the time 
being, let us assume that this matrix has been constructed. The usual prescription 
proposed in [2] to obtain the Green's function is generalized in the present setup as 
follows. Rewrite the action (2.5) by freeing the fields from derivatives. After 
inserting (2.9a) the action can be written as follows 

c^jE.O.M.]^ + ^[2Ajj^_ k ^ k J + B\j^(] 
= J dk> <pL k ru(k,z)<p J k \* , (2.11) 

6 K.L. would like to thank G. Moore for a discussion on this point. 
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where Zh and = stand for the limiting values of z at the horizon and the boundary 
respectively. In the last line we passed to the on-shell action and defined the flux 
matrix 

F(k, z) = 2F ] A H F' + F ] B ] F . (2.12) 

From here, the natural generalization of the original Minkowskian AdS/CFT 
prescription [2] is just 7 

G iA k ) = - lim Fu(k, z A ) ■ (2.13) 

Strictly speaking we have derived this relation for "positive" k = fc>. However (2.13) 
extends smoothly over to "negative" k = fc<. To see this, one has to start however 
from the same bulk action (2.5) and instead free from derivatives, hence making 
use of the appropriate equations of motion. The boundary action then adopts exactly 
the same form as in (2.11) with the replacement k — > —k in the integrand. Given 
the conjugation properties of the matrices A, B and F under change of sign in k this 
is consistent with the required property of retarded Green's function (see appendix 
A) 

GfA-k) = Gfj(ky. (2.14) 

To conclude this section, let us mention that the definition of the Green's function 
as given by equation (2.13) is still somewhat incomplete. The bilinear action we 
wrote will generally present divergences at the boundary that must be regularized 
by adding appropriate covariant counterterms to the action. These counterterms 
change the definition of the flux matrix by 

Fu{k, z k ) -> Tij{k, za) - Fct,u{k, z K ) , (2.15) 

which in the limit za — » gives a finite answer. The exact form of the terms to be 
added depends on the theory under consideration, and for the case studied in this 
paper, the appropriate expression for T c t,ij{k, Za) can be found in section 3.5. 

2.3 The Green's functions as bulk Noether currents 

Due to the arbitrariness of the (p^, equation (2.6) implies the following 



[-2(A H FJ + 2B A F' + (2C H - B V )F] U = . (2.16) 

We can multiply this from the left by F J M{—k, z) = F M J (k,z)* = Ft M J , thus 
obtaining the following matrix statement 

-2F\A H F')' + 2F ] B A F' + F\2C H - B ] ')F = . (2.17) 



7 In [16] the same authors also deal with a mixed operator situation. We seem to disagree with 
their prescription (see eq. (4.26) there) in which diagonal and off diagonal components are treated 
on different footings. 
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Now proceeding as before, obtaining the equations of motion by varying inserting 
(2.9b), and then contracting from the right with F we end up with the adjoint version 
of (2.17) 

-2(F n A H )'F - 2F' ] B A F + F ] (2C H - B')F = . (2.18) 
Subtracting (2.18) from (2.17) and using (2.12) we obtain 

±(F(k,z)-^(k,z)) LJ = 0. (2.19) 

The fundamental reason for the conservation of so many quantities is because we 
work at the bilinear level of the action. Once we have written this in terms of the 
and <3>_fc we can assume the positive and negative frequency fields to be independent. 
The complete bulk action (2.5) can be written as 



S = I dk> 



Y>LkS?M\ > (2-20) 
where Sjj is the matrix of 1-dimensional action functionals 
S$ = J dz Cf)(z) = J dz [2F'*A H F' + F^BF' + F'^B^F + 2F*C H F] IJ , (2.21) 

for the N 2 one-dimensional "operator mixing fields" F 1 j(k, z). For each /, J we find 
a U(l) symmetry F 1 j(k, z) — > e iaiJ F 1 j(k, z). Hence, for each k, we obtain a matrix 
of Noether currents 

/(*) (A - ( ,-, dC^ MN (z) j dC^ MN {z) fJ 

j mn\ z ) - (+V — q^Tj — * J + l-V dF>]Ij h J 

= i(2F'^A H F + F ] BF) MN - i{2F^A H F' + F ] B ] F) MN 

= -i{F-J*) M N. (2.22) 

The evaluation of this current in the case under consideration in this paper can be 
found in section 3.4. 

Notice that in terms of Green's functions the ^-independent quantity is precisely 
the matrix spectral function, p{k) = i{G R {k) — G A {k)) (since G A {k) = G R {k)\ see 
appendix A), which therefore turns out to be an RG flow invariant quantity. In fact, 
the evaluation of the analog to this current for different systems was the tool used to 
study the phenomena such as graviton absorption prior to the celebrated prescription 
of Son and Starinets [2] (see for example [17, 18]). 

2.4 Quasinormal modes 

In gravitational scenarios there are different ways of defining quasinormal modes, 
depending on the boundary condition we impose on the fluctuations of the fields. 
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This degeneracy of boundary conditions is lifted in AdS/CFT by stating that the 
quasinormal modes relevant for the holographic interpretation must correspond to 
poles of the holographic Green's functions [19]. Indeed, the Green's function matrix 
we have just defined will in general be a meromophic function of frequency and 
momentum. We therefore define the quasinormal modes of the coupled system as 

quasinormal modes <-* poles of G R (k) . (2.23) 



In the analysis below we will focus on real positive values of the spatial momen- 
tum, so the quasinormal modes will be given by complex frequencies uj n = Q n + iT n , 
where — T n gives the damping factor of each mode. The presence of modes with pos- 
itive T n signals instabilities of the system, as for these modes their amplitudes grow 
with time. These unstable modes define tachyonic instabilities. Once the quasinor- 
mal modes are known one can express the meromorphic Green's function as a sum 
over poles plus an analytic part 8 

G>, q) = £ nniq) ( ; + T(u, q) , (2.24) 

n=l v ' 

where TZ n (q) and T(ou, q) are NxN matrices with analytic components. As explained 
in appendix A the spectral function and the causal Green's function are related to 
one another as follows 

p(u, q) = z [G R (u, q) - G R (u, q) f ] . (2.25) 

Provided we are exploring the real u axis, the full spectral function matrix can be 
expressed in the Breit-Wigner form as 

^ (u - n n )( Im [K n ] + Im [TZ n } T ) + T n (Re [K n } + Re [JZ n ] T ) 

Re [p] = ~ ^ (u - n n y + vi ' (2 ' 26a) 

j r I V- (w - n n )(Re [n n ] - Re [TZ n } T ) + T n { Im [K n ] - Im [TZ n } T ) 



n=l 



plus the contribution from the analytic part. Notice that, consistently, the diagonal 
terms of the spectral function are real. 

From the discussion above we see that to determine the properties of the Green's 
function it is important not only to find the position of the quasinormal modes, but 
also their residues. In appendix B a numerical recipe to obtain the residues can be 
found. We have checked that the position and residues of the quasinormal modes 



8 We assume here that simple poles are the only non-analyticities in the holographic Green's 
functions. Although there is no proof so far in the literature this seems to be the case for the 
non-extremal asymptotically AdS black holes of relevance here. 
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obtained by the methods given here and by fitting the spectral function in the real- 
frequency axis to a Breit-Wigner function give compatible answers, up to a smooth 
analytic function of the frequency. 

When dealing with a parity-invariant system, the retarded Green's function sat- 
isfies 9 GfAuj, q) = o-j(jjG^(— lj, q)* implying that poles must come in pairs such that 
there is a relation between them 

^m,ij(q) = -^jK,ii(q) , w m (q) = -u£(q) , (2.27) 

for fixed m and n. This relation classifies the quasinormal modes into two different 
types. On one hand, when n ^ m we observe that each mode has a "dual" mode 
with position and residue given by the former relation. Obviously, there are an even 
number of these, half with positive real part and half with negative real part. They 
are responsible for the quasiparticles observed in the spectral function, as will be 
clarified later. On the other hand, when n = m we find purely imaginary modes 
with the corresponding residue matrix being purely imaginary. Some of the modes 
may satisfy the limit 

lima> n (q) -> . 

q^O 

These are the only modes that survive at long wavelenghts and long times, therefore 
we call them "hydrodynamic modes". From them one can extract all the relevant 
information about the hydrodynamic properties of the system. It may be that two 
modes of one of the classes stated can recombine, becoming two modes of the other 
class. When this happens, both modes must have zero real part and their residues 
must also be purely imaginary. 

Former studies of the residues [12, 20] by fitting the spectral function resonances 
to a Breit-Wigner function did not take into consideration the possibility of a complex 
residue. The complex residue acts by introducing a phase in the quasinormal mode 
that shifts the position of the maximum of the spectral function in the real cu-axis 
with respect to the position Q n of the pole. 

2.5 Adapting the prescription to numerical solutions 

Except for some simple cases one does not expect to find an analytic solution to the 
iV coupled equations (2.6), and therefore it is not possible to extract the solution 
matrix F(k, z) analytically. It follows that we are forced to give a prescription to 
calculate this matrix from numerical results. At the level of fluctuations, we work 
with a bilinear action, and hence equations of motion are linear and second order. 
Hence, on general grounds, we expect to find a basis of 2N solutions. To obtain 
any of these solutions we must supply boundary data at a given point from which 

9 <x; = ±1 is the charge under parity reversal of the operator Oi. See appendix A for details. 
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integration starts. Whenever a black hole is present in the geometry, the event 
horizon is the convenient position at which to impose boundary conditions. This is 
so because automatically we can halve the number of basis solutions by demanding 
"in-going" boundary conditions at this point, which is what leads ultimately to the 
computation of a retarded Green's function [2]. 

Having fixed N boundary conditions, the other N correspond to normalizations 
of the fields. We can select among N independent iV-tuples that can be chosen to 
be 

*(.) = (* " Zh)~* (e{ a) + 0{z - z H )) , (2.28) 

where ro oc uj/T is a dimensionless frequency, weighted by the Hawking temperature 
of the black hole 10 . The N linearily independent vectors ei a \ can be chosen to be 

efo = (I,-- - ,1), (2.29a) 

e U = (V" >s^L/--- « = {2,--- ,N}. (2.29b) 

a th 

Therefore, we have given N sets of independent boundary conditions at the horizon. 
We can perform a numerical integration for each set and obtain N independent 
solutions that extend in the range z G (z A , z^). Let us call these, the IR-normalized 
solutions, and arrange them in a matrix, H(k, z), in such a way that the J th solution 
(<3>(j), ^(j)) appears as the J th column, i.e. 

H I j(k,z) = $ I {J) (k,z). (2.30) 

Any "in-going" solution, can be written as a linear combination of these N indepen- 
dent solutions. In particular the matrix F{k, z) of UV -normalized solutions must be 
linearily related to H(k, z). Since at the UV cutoff, by definition, F{k,Z\) = 1, the 
linear relation must be 11 

F(k,z) = H(k,z) ■ Hik^ZA)- 1 . (2.31) 

In general we will take the limit z\ — >• to evaluate the expressions at the boundary. 
As stated before, the Green's function is given by (up to regularizing counterterms) 

G R {k) = - lim F(k,z A ) = - lim (2A H (k, z A )F'(k, z A ) + B\k, z A )) , (2.32) 

10 For example, in the quenched D3/D7 system we are going to study in later sections, ro = ■ 
n The case of N uncoupled fields is automatically included. In this case, the matrix F is by 
construction diagonal for all z, and takes values 

F I J = di a g[$ 1 (z)/$\z A ),--- ,$ N (z)/$ N (z A )] . 



- 11 - 



where we have taken into account the UV normalization of the matrix F(k, z\). Now, 
after having made sure that the behaviour close to the boundary (z\ — ► 0) is 

H T j(k, 0) ~ A(k) T j + z A +~ A -B(kY j + ... , (2.33) 

with A(k) and B(k) the connection coefficient matrices, we can insert this into (2.31) 
and (2.32) and get 12 



J 



G R {k) T j = - lim 2(A'-Ai)z* + A ~ 1 fA H (fc, ^(fc^fe)- 1 ) + flt(Jfc, z A ) 

(2.34) 

Note that the non-analytic (in k) behaviour comes from the Au terms in the action. 
The Bjj terms will give analytic contributions to the Green's function. 

Moreover, notice that G R (k) is ill-defined whenever detA(k) = 0. From equa- 
tion (2.33) we see that the Green's function has poles whenever the inverse matrix 
H(k,z\)~ l does not exist, which is consistent with the discussion in section 2.4 by 
equation (2.31). Under the present construction, this is equivalent to demanding 
that the determinant of H vanishes at the cutoff 

det[H(k n ,z A )] = 0, (2.35) 

which is a very convenient operational statement for determining the position of the 
quasinormal modes in the complex frequency plane numerically. With it, one can 
track the position of quasinormal modes whose effect cannot be observed (or even 
guessed) in the spectral function due to their being too far down into the imaginary-^ 
axis, or the associated residue's value being small. 



3. Example application: D3/D7 probe fluctuations at finite 
baryon density 

For completeness, in this section we describe a system consisting of a set of Nf 
probe D7-branes in the background of a stack of N c non-extremal D3-branes with 
Nf N c . Our notations and conventions will be as in [13]. We will apply the 
formalism developed in the previous section to compute the quasinormal modes and 
spectral functions in section 4. 

3.1 Background 

In the framework of the AdS/CFT correspondence, the retarded correlators we are 
interested in, G R (k), can be obtained from the perturbations of a U(l) gauge field 



12 The reader will recognize here the generalization of the G R ~ B/A rule of thumb put forward 
in [21]. 



dual to the electromagnetic current on the boundary. The relevant holographic de- 
scription is provided by an AdS geometry with a non-extremal horizon and embed- 
ded probe branes. The baryonic U(l) symmetry is the abelian center of the natural 
U (Nf) global symmetry present on a stack of Nf coincident D7-branes. For the case 
of interest here, namely D3/D7 configurations, the dynamics of this gauge field is 
encoded in the action for the probe Dg-brane 

S = -N f T D7 J d 8 £ v/-det(# + 27ra'F) + S wz . (3.1) 

The second term on the r.h.s. stands for the Wess-Zumino term which will not 
contribute to the equations of motion for the background or the fluctuations. Td 7 = 
l/((27rl s ) 7 g s l s ) is the D7-brane tension, g s is the string coupling constant and g is 
the pullback metric induced by the relevant background. As for the background, we 
will be dealing with the near horizon limit of a stack of non-extremal D3-branes 

ds 2 = H-^i-fdt 2 + d4) + H 1 ' 2 + P 2 dn 2 ^ , (3.2a) 

where x 3 = (x l ,x 2 ) x' i ) and 

H W={j) ; M = i-(ff- (3-3) 



L = 4:ng s N c lg and the Hawking temperature is given by 

(?) ■ <-> 

The D3/D7 intersection is summarized in the following array 



T = — f p ° S 
ttL 



0123456789 
D3 : x x x x 
D7 :xxxxxxxx 

where the probe D7-branes wrap a 3— sphere in the directions transverse to the 
_D3-branes, so it is convenient to write the metric of S 5 in adapted coordinates, 

dQl = d6 2 + sin 2 6 dQ 2 3 + cos 2 6 d<f) 2 . (3.5) 

Setting ip = cos 6 the classical Dg-brane embedding may be specified by a dependence 
if} = if>{p)- For numerical analysis we have changed to a new dimensionless radial 
coordinate u related to p by 

»=f?V. (3.6) 



-13- 



In terms of this, f(u) = 1 — u 2 , the horizon (boundary) lies at u — 1 {u — 0) and 
H(u) = j^fjjyi- Specifiying the _D7-brane embedding through ip = ip{u) the induced 
metric takes the form 

>9 (vrTL) 2 , - , 9N L 2 (l-ib 2 + Au 2 f^' 2 ) , 9 9N „ , s 

= ^7^(-/ rft + rfx s) + 4 ^j (1 _^ 2) du + L (! - ^ 2 )^3 • (3-7) 

The generalization of the previous setup for finite baryon density was investigated in 
[9, 10]. The relevant bulk degree of freedom dual to the baryon chemical potential is 
the Aq component of a U(l) gauge field on the world-volume of the D7-brane. The 
background profiles for ip(u) and Aq(u) are obtained by solving the Euler-Lagrange 
equations of the Born-Infeld lagrangian in (3.1). The gauge field A (u) obeys a 
conservation equation owing to the fact that it enters the action purely through its 
derivatives, so its solution can be expressed in terms of a constant of integration d 
as follows 



4(«) = -d— V (3.8) 
4a ^ 2 {^ + d 2 vfi) 

where tjj = ^/l — ip 2 . This constant of integration is related to the electric displace- 
ment n q by 

n q = (3.9) 
8a' 

If one wants to integrate this expression further, the correct condition to impose 
by regularity is that the gauge field at the horizon vanishes. With this regularity 
condition the quantity fi = Ao(0) is holographically identified with the chemical 
potential [22] . We see that in the limit of vanishing baryon density d — > we obtain 
vanishing chemical potential (although there is a region of the phase diagram for 
which this does not hold for sufficiently large quark mass [10, 23]). 

The equation for ip{u) gives 

/ AfftiP'Jfo + dW \ *l> (3^ 4 + 4uW'W - d 2 u 3 )) 
d u [ —. V ^= +^ _ J =0, (3.10) 

\uJ^{i) 2 + Au 2 fip' 2 ) J u 3 ^ 6 (^ 6 + d 2 u 3 )(4> 2 + Au 2 f*l)' 2 ) 

and cannot be solved analytically. Close to the boundary this equation reads d u (Aip'/ u) 
—3ip/u 3 and its solution behaves as 

M U ) ~ ^-u 1 ' 2 + -^u 3 / 2 + 0(u 5 / 2 ) , (3.11) 
whereas close to the horizon 13 

* ^ ~ t ( M \~f 0) ]~^ " «) + 0(1 - uf . (3.12) 
o (1 - ipffi + d 2 

13 In the presence of a non-zero baryon density all the embedding profiles of the D7-brancs reach 
the horizon [10]. 
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The series near the horizon depends only on one parameter ip e [0, 1) and can 
be used to integrate numerically towards the boundary. Once this integration is done 
we can read off the asymptotic behaviour and extract the boundary quantities m(ipo) 
and c{ipo) from equation (3.11). These constants m and c parametrize respectively 
the quark mass and something we loosely refer to as the quark condensate [6, 24, 25, 
26, 27], 

M q = -\fXTm, (O) = -\^\NfN c T z c, (3.13) 
2 8 

with A = gy M N c = 2iig s N c , the 't Hooft coupling. The operator O is a supersym- 
metric version of the quark bilinear 

O = + $ f X$ + M q &$ , (3.14) 

with X one of the adjoint scalars. A precise definition can be found in [10]. The 
3-area of the induced horizon (per unit 3-dimensional Minkowski space volume) is 
controlled by ip Q 

A H = 27r 2 (7rTL 2 ) 3 (l - ^Plf' 2 . (3.15) 

We expect this quantity to govern the rough shape of the peaks of the spectral 
function with larger widths for larger induced horizons. We are going to refer to this 
quantity several times in section 4, so we find it convenient to plot it in figure 1. 

3.2 Fluctuations 

We will consider perturbations of the world-volume fields that depend only on the 
RG flow coordinate u and the Minkowski coordinates a; , a: 1 , thus not considering 
any dependence on the internal coordinates wrapping the S 3 . 

i/>(u,x) -> + ee-^^'^H, (3.16a) 

A^u,x) - A» +ee-^°-^2l,(u) . (3.16b) 

With this we can expand the DBI lagrangian in powers of e 

C = C + eC 1 + e 2 C 2 + --- . (3.17) 

Upon imposing the equations of motion for the background fields, L\ vanishes and 
the linearized equations for the perturbations can be derived from the quadratic 
piece. 

The fields 2to, 2li and ^ form a coupled system of differential equations. At first 
sight this seems somewhat surprising since the scalar field is uncharged under the 
corresponding U(l) gauge symmetry. We must remember however that in the case of 
a non-trivial gauge field background we are really dealing with charged flavor probe 
branes and that \l/ parametrizes the deformations of the flavor branes around its 
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Figure 1: Normalized induced horizon area on the D7-branes as a function of the quark 
mass and the baryon density. Red lines mark regions of equal induced area with the corre- 
sponding value of ipo specified. For d < 0.00315 the curve is multivalued close to m = 1.3, 
signaling a first order phase transition. Furthermore, it was recently shown in [30] that 
an unstable quasinormal mode with positive imaginary part of the frequency exists in that 
region. We will however not consider it in the current paper. 

equilibrium configuration. Thus, if we deform the probe branes the charge distribu- 
tion on them will also experience induced forces since now it is not in equilibrium. 
Therefore, scalar field fluctuations, necessarily will also induce fluctuations in the 
charge density. We might think of the the scalar field as carrying multipole charges 
with respect to the gauge field. Indeed, upon expansion of the DBI action, couplings 
of the scalar field to the field strength tensor do appear. They are caused by multipole 
moments of the charge distribution and vanish therefore at zero momentum. Since 
the scalar does not carry monopole charge it is still possible to rewrite the equations 
of motion using a single propagating gauge invariant combination El where 



is the longitudinal (parallel to the fluctuation) electric field. There are no couplings to 
the transverse fields Et = cu2l2,3- The equations of motion can be found in appendix 



The definition of quasinormal modes as the zeroes of a determinant spanned by 
the field values at the boundary has been used before in [14, 28]. As explained there, 
it is also possible to avoid the introduction of gauge invariant fields and work directly 



E L = q% + uQL 




D. 
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with the gauge fields. In such a case pure gauge configurations need to be taken into 
account in order to obtain a maximal set of linearly independent solutions. In the 
case at hand it is however easier to work with the gauge invariant electric field (3.18). 

In order to solve the coupled equations of motion we must setup the boundary 
conditions for the coupled system. The lore is that for retarded Green's functions 
we must select incoming wave boundary conditions on the black hole horizon. The 
generalization of the usual Frobenius expansion near the horizon is straightforward at 
first order, and in our case gives the usual regularizing factors $(w) — > f(u)~ l ^r&(u). 
Now the regular coupled system of equations can be numerically integrated from the 
horizon towards the boundary to obtain the Fourier bulk modes *&k( u ) an d El^{v) 
with k = (u, q, 0, 0). 



3.3 Green's function from fluctuations 

As mentioned before, in principle we are perturbing the gauge fields 2t M and the 
scalar field However, gauge symmetry and the fact that we only have rotational 
invariance in the thermal vacuum implies that the relevant fields are the gauge in- 
variant combinations El = g2lo + cj21i, = UJ^U (i = 2,3) and At this stage we 
must write the boundary action in terms of the ^fe(w) and E L k (u), E Ttk (u) degrees 
of freedom. The neatest strategy is to write the bulk bilinear action in the form 
given in equation (2.11). Hence we must proceed by writing the Fourier transformed 
action in the form (2.7) and extracting from it the explicit values of Ajj and Bjj. 
Details can be found in appendix C. 

From the form of the boundary action we expect a structure of retarded corre- 
lators given in terms of these fields as follows 



G R (E T ,E L ^) 



( (E T E T ) \ 

(E L E L ) (E L $) 
\ (ifE L ) (M) J 



(3.19) 



From this matrix, it is straightforward to obtain the correlators for gauge fields. In 
fact, defining the polarizations 



G R {E T ,E L ,V) = 



/ n r (fc) 

LU 2 




V 





U L (k) 



2 — q 2 
(k) 



W — q 





n L *(fc) 



\ 



(k) 



(3.20) 



/ 



all the relevant information is contained in this set of functions. At q = rotational 
invariance is restored, implying U L (u>,0) = n T (u;,0) and n L *(u;, 0) = 0. Also from 
the requirement that the Green's function is regular on the light cone we must find 
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for k 2 = that U L (k) = H Li ^(k) = 0. The usual correlators for conserved currents 
are obtained from here through the introduction of the relevant kinematical factors, 
which are found by using the chain rule 



5W 



(3.21] 



where i,j = T,L and the transverse and the longitudinal projectors are defined in 
the standard way (see [15] for example) For /c M = (u, q, 0, 0) this leads to the only 
non-vanishing components 



C^=IT>,g), * = {2,3} 



(3.22) 



2121 



I 2 ttLs . „a ^2121 ~1 U} t-,L 



lu 2 — q 2 



n>,?), c£f 



n>,?), cS 



co> 2 — g 2 



or 



lu 2 — q 2 



U L (u;,q) 



and for the Green's function 



521" 



again with i = T, L we obtain from (3.19) and (3.20) 



21* 



=n 



21* 



v 7 ^ 



=n 



L* 



(3.23) 



3.4 Conserved current 



For the D3/D7 system we can evaluate the Noether current at the horizon, reading 
the matrix expressions given in appendix C. We can use the IR- normalized matrix of 
solutions, H(k,u), to perform the derivatives and then evaluate them at the horizon. 
The holographic information of the system enters through the factors of H _1 (k, 0) in 
the definition of the UV-normalized matrix of solutions, F(k, u), which are the ones 
entering naturally in the definition (2.22). 

The A H (k, u) matrix can be shown to behave near the horizon as 0(l — u) in the 
diagonal terms and 0(1 — u) 2 in the off-diagonal ones. The B^(k, u) matrix has a null 
diagonal and the off-diagonal terms behave like 0(1 — u). Therefore, the evaluation 
of the matrix of Noether currents gives 



J(k) 



lim 

U->1 



/ 



(2tiTuj) odc F(k, u) 



(2-kTuj) 





(txTL 



2\4 



F(k,u) 



(3.24) 



where udc = _ 4>o) 3 + <^ 2 * s the DC conductivity of the system. 
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As stated above, this quantity is identified with the spectral function of the 
system. It is straightforward to show from the former expression that the diagonal 
entries of uj pijiui) are positive. 

Numerically we have checked that the antihermitian part of the flux J-(k, z) 
is independent of the radial variable in the parameter region where numerics are 
to be trusted, in full agreement with equation (2.19). We also checked that the 4 
independent components of this antihermitian matrix are given by expression (3.24). 

3.5 Regularized action 

The counterterms needed to regularize the D3/D7 quenched system were obtained 
in [29] and can be expressed as 

J d 4 xS ct = - J d 4 x^^(l - ^ 2 ) 2 , (3.25) 

where is the euclidianized boundary metric and i/j is the embedding profile. When 
one perturbs this profile by considering 

4){u) -> rj)(u) +ee -iK-^)^$( u ) , (3.26) 

where the normalization factor yfu has been taken into consideration (see appendix 
C), then the counterterm action can be expanded in powers of e, which effectively 
marks the number of perturbation fields 14 . At second order the counterterm enters 
the definition of our boundary action, which is now defined as 

S = Jd~k> (2Afj¥_ k &i + B\j¥_M - 2S l!: A> ,vI/„) , (3.27) 

where ^>k = &{ =2 goes to a constant at the boundary. It is easy to check from 
equations (C.3) that the A H matrix is regular at the boundary, whereas the B< 
matrix reads 

(0 \ 

B\u -> 0) = I (ttTL 2 ) 4 + 0(1) (3.28) 
Close to the boundary the counterterm quadratic in the fluctuations gives 

Sat^u - 0) = + 0(i) (3.29) 

so the contribution to the boundary action is 

'o 

£ f - ^ - 2 I (ttTL 2 ) 4 I + Oil) = Bl gular , (3.30) 







2u 



14 In this subsection we recover the barred notation of section 2.1. 



-19- 



and the Green's function is divergence-free with the usual counterterms. It is worth 
noting that the added counterterms affect only the hermitian part of the flux matrix, 
because they enter through the real part of a diagonal component, this means that the 
spectral function for this system is insensitive to the presence of these counterterms, 
as is the position of the quasinormal modes. This is consistent with the fact that the 
spectral function is ^-independent. 



4. Analysis and discussion 
4.1 The mixing mechanism 

There are three independent limits in which one can find that the system considered 
in the previous section decouples, these are the massless quark limit m — > 0, the null 
momentum 15 limit q — > and the zero baryon density limit cZ — ^ 0. When none of 
these limits is taken we have to face the presence of coupled fields. A question arises 
naturally, how does the mixing appear from the point of view of the quasinormal 
modes? 

One convenient way to find an answer to this question is consider first the decou- 
pled case. Taking one of the decoupling limits one can study either the longitudinal 
electric field sector or the scalar sector without considering the other. Then, numer- 
ically one can find the quasinormal modes. This was done in [30]. It was observed 
that the quasinormal modes of both channels do not coincide at finite temperature. 

Returning to the coupled case (for example, evolving the parameters slowly from 
a decoupling limit), the system becomes coupled and we cannot talk anymore about 
poles associated to the longitudinal vector sector or to the scalar sector. The poles 
are collective properties of the system. Despite this, one would like to understand 
how these collective modes can be categorized in the decoupling limits and identified 
with one of the two channels under consideration. 

In figure 2 we compare the finite temperature contribution of the longitudinal 
electric field component of the spectral function 

NfN ° T2 -PaM = i [G R (tt>) - G R \rv)} 1 x - I ^^2n 6(ra 2 - q 2 ) , (4.1) 



4(ro 2 - q 2 ) 

with the position of the quasinormal modes of the system for the same parameter 
values. 

We see that when we deal with the coupled system the poles appear in prox- 
imate pairs with similar values Q n . One can calculate for each of these modes its 
corresponding matrix of residues. Taking one of the parameters (m, q or d) to evolve 
towards the decoupling limit and studying how these matrices of residues change, we 

15 Throughout this section we will work with the dimensionless momentum (to, q) = £'t ■ 
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PaM // Im[l»] 
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(b) 

Figure 2: Example of the position of the quasinormal modes with positive real part (red 
points, scale in negative axis) and the corresponding finite temperature contribution to the 
component of the spectral function associated with the longitudinal electric field propagator 
(continuous line, scale in positive axis). Notice that in figure (a) only half of the poles seem 
to contribute to the spectral function, this is because the other half have a small residue. 
In figure (b) we plot a detailed version of a spectral function where all the poles have an 
observable contribution. These plots are for m = 0.01, d = 2 and q = 0.2 and q = 2.2 
respectively. 



observe that in the case where the parameter is small (this is, when the system is 
weakly coupled) the matrices of residues for the proximate pairs of modes tend to 

which is exactly what one expects to find if the system were decoupled. 
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This implies that in the decoupling limit we can state that the quasinormal 
modes decouple by means of the matrix of residues, and one can then associate each 
of these poles either to the longitudinal vector channel or the scalar channel. With 
this interpretation we are able to recover former results found in the cited literature. 
When we are close in the parameter space to the region where the coupling is small, 
the shape of the spectral function resembles closely the spectral function of the 
decoupled cases (see figure 2(a)). However, once we further probe the parameter 
space, the magnitudes of the residues associated to the proximate pairs of poles 
become similar to each other, and thus the peaks of the spectral functions contain 
a more complicated structure, i.e., each peak in the spectral function has a deeper 
structure given by the contribution of two poles 16 (see figure 2(b)), one of which can 
be linked to the scalar channel in a decoupling limit, and the other to the longitudinal 
vector channel. 

Yet another way to see how the mixing appears in the system is to focus on 
expression (3.24) for the spectral function. In a decoupling limit the matrix F(k,u) 
is diagonal, thus giving a diagonal spectral function, each term of the diagonal corre- 
sponding to each of the uncoupled channels. Correspondingly, we have two indepen- 
dent Green's functions. The matrix F(k, u) is sensitive to the bulk of the holographic 
geometry, and when the system departs from the uncoupled case, this matrix will 
notice the mixing in the bulk of the two fields, and will no longer be diagonal. This 
means that the spectral function is now given by a 2 x 2 matrix, and the same holds 
for the Green's function of the system. 

4.2 Field theory interpretation 

In order to elucidate the connection between field theory effects and holographic 
mixing, we discuss the significance of mixing fields in the bulk from the boundary 
field theory point of view. 

Renormalization Field coupling in the bulk means that a single field <&{(z) sources 
a linear combination of all operators at the cut-off z\. Any field may be expressed 
through the bulk to boundary propagator Ff and the boundary data (p{ as seen from 
equation (2.9a). Hence the bulk to boundary propagator Ff describes the behavior 
of an operator Oi under the RG-flow along the radial coordinate z, i.e. the operator 
renormalization Q™ nmmahzcd — Z h a O\ &xc . This is in analogy to the renormalization of 
fields in ordinary quantum field theory, where we have ^ renormallzed = Zijj ha,Te . 

Furthermore, the coupling of the gravity fields introduces renormalization 
corrections from all operators to the two-point functions. That means that a single 
operator Green's function ([Oi, Oj\) with a fixed / is in general renormalized through 

16 In reality each peak has contributions coming from all the quasinormal modes, but these con- 
tributions die away as (lo — f2„) _2 . 
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contributions of all operators. This is in analogy to the loop corrections describing 
renormalization in ordinary interacting quantum field theory. Inside these diagrams 
all particles (with appropriate interaction vertices) of the theory may appear, just like 
in our case all the operators may appear. However, in contrast to ordinary quantum 
field theory here we do not have to compute higher loop orders in order to get the 
full action of the renormalization group. Due to the gauge/gravity correspondence 
the exactly renormalized field theory result is (in our semi-classical approximation) 
already encoded in the leading order gravity solutions. So this effect is an exam- 
ple of how quantum effects such as renormalization of a quantum field theory are 
holographically encoded in the dynamics of a purely classical bulk theory 17 . 

4.3 Hydrodynamic regime 

The vector field on the brane corresponds to a global U(l) symmetry in the dual field 
theory. At finite temperature the global symmetry has to give rise to a hydrody- 
namic mode since a conserved charge can not be dissipated away but diffuses slowly 
through the medium. In other words we expect to find a quasinormal mode with 
a hydrodynamic dispersion relation such that \imu(q) = 0. Furthermore for small 

momenta the dispersion relation has to take the form of a diffusion kernel uo = —iDq 2 
where the diffusion constant now depends on the ratio of quark mass to temperature 
and the baryon density d. 

Upon increasing momentum, higher powers appear in the disperion relation giv- 
ing rise to higher order hydrodynamics. If we increase the momentum still further, 
we expect however a crossover from diffusive regime to a reactive regime 18 . More 
precisely we expect the hydrodynamic diffusion to show up as a purely exponential 
decay in time where as at smaller wavelengths we expect to find a slowly decaying 
oscillating behavior. In the holographic context this crossover has first been dis- 
cussed in [32] by studying spectral functions. This crossover can however also and 
more directly be addressed in terms of the quasinormal modes. It has been observed 
in [33, 34] that for the longitudinal R-charge current and the shear channel in the 
AdSs black hole background, that there exists a certain critical momentum value 
from which it is not anymore the purely imaginary diffusion mode that dominates 
the long time behaviour of the system. More precisely, from that value of the mo- 
mentum on, the imaginary part of the first non-gapped quasinormal mode is closer 

17 Note that in the case of the dilepton production rates only the diagonal current-current corre- 
lator contributes, since it is only the current that couples to the intermediate off-shell photon that 
decays into the dilepton in the final state [31]. This is true even in our case where the longitu- 
dinal current components mix with the scalar. Nevertheless, the scalar contributes virtually since 
the scalar fluctuations are necessarily switched on in the bulk and the scalar quasinormal modes 
influence the shape of the current spectral function if the mixing is large enough. 

18 In the literature this is sometimes also called a hydrodynamics to " collisionlcss" or "quasipar- 
ticle" crossover. 
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to the real axes than the diffusive mode. In terms of time development this shows 
up as a change from a purely damped decay to a slowly decaying oscillation since 
the first gapped mode also has a non-vanishing real part of its frequency. In [30] it 
was shown that this behaviour also holds for the longitudinal sector of the vector 
field on the D7-brane at finite quark mass. Moreover, the change in behaviour in the 
time domain has been explicitely verified recently in [35]. We can therefore identify 
the value of the momentum where the purely imaginary diffusion mode crosses the 
imaginary part of the non-hydrodynamic mode, moving from the diffusive to the 
reactive regime. 

We would like to investigate how this crossover takes place in the case with non- 
vanishing gauge field on the D7-brane. Naively one might expect that nothing new 
would happen compared to the case without baryon charge. However, now we have 
to take into account the mixing of the longitudinal vector channel with the scalar one. 
An important feature of the scalar sector quasinormal mode spectrum is the existence 
of purely imaginary poles, as has been shown in [30]. For vanishing or small densities 
these modes are responsible for the appearance of a tachyonic instability at low 
temperature to mass ratios. There exists a rather small critical density d = 0.00315 
above which the system becomes stable. If we switch on momentum we know that 
the scalar channel and the longitudinal vector channel mix. If there are now two 
neighbouring purely imaginary quasinormal modes in the spectrum, it may happen 
that they combine and move off the imaginary axis developing non-vanishing real 
parts. In fact, this is the way the crossover from the hydrodynamic regime happens 
in AdS4 [36] and on probe D5-branes corresponding to defects in the CFT [20] . 

In figure 3 we have plotted the dispersion relation for the relevant modes at a 
fixed quark mass m = 1 but for different baryon densities. The left plot (a) shows 
the imaginary parts of two modes for a rather low density d = 0.01. There is a purely 
imaginary, hydrodynamic mode representing the diffusive behaviour of the baryon 
charge. The other mode is gapped and has a real as well as an imaginary part. Only 
the imaginary part is shown since this determines the decay time. As we can see, 
there is a critical value of the momentum at around q = 0.8257 where the diffusive 
mode crosses the non-hydrodynamic mode. From that value on the response of the 
system is dominated by the first non-hydrodynamic mode and we may say that it 
has entered the reactive regime. No other mode is visible at this small density in this 
region of frequency and momentum space. This is qualitatively the same behaviour 
as in the zero density case [30]. At larger momenta the purely imaginary diffusion 
mode might pair up with another purely imaginary mode. For the crossover to the 
quasiparticle regime this is however not relevant, since the long time behavior is 
already dominated by the lowest quasinormal mode shown in 3(a). 

The middle figure 3(b) shows the situation at a higher density of d = 0.545. 
Now we see three modes. There is the diffusive mode, the first gapped quasinormal 
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(a) 



(b) 




(c) 



Figure 3: Crossover from the diffusive to the reactive regime in terms of quasinormal 
modes. Two different mechanisms of how this crossover happens can be seen. At small den- 
sity (a) the hydrodynamic mode crosses the imaginary part of the lowest non-hydrodynamic 
mode. At large density (c) the hydrodynamic mode pairs up with another purely imaginary 
mode and moves off the imaginary axes as a pair with non-vanishing real frequencies. In 
between (b) the three imaginary parts of the modes meet at a single value of momentum. 



mode that also has a real part (plotted as a dashed line) and there is now a second 
purely imaginary mode. As we increase the momentum all three lines meet in a 
single point at q = 0.73 and for larger momentum only two lines are visible. This 
corresponds to the fact that the two purely imaginary modes have combined into 
a pair of quasinormal modes with non-vanishing real parts. We can identify the 
momentum where all three lines meet at the point where the crossover from the 
hydrodynamic to the quasiparticle regime takes place. 

The lower figure 3(c) shows the same modes now for a rather high density d = 2. 
Now the two purely imaginary modes combine first into a pair of quasinormal modes 
with non-zero real part . We can identify the crossover now with the momentum 
at which this pairing of the purely imaginary modes takes place. In figure 3(c) this 
takes place at q = 0.35. 

From the the dispersion relation we can also compute numerically the diffusion 
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constant D as a function of density and temperature. This problem has been ad- 
dressed before and in [37] a general formula in terms of the background fields has been 
derived. Our numerical results for the diffusion constant are in very good agreement 
with this formula 



D- aDC 
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(4.3) 



with odc the DC conductivity, \ the susceptibility and A, S and 7 defined in ap- 
pendix D. 

4.4 Quasiparticle regime 

In this section we will focus on the regime where peaks on the spectral function 
can be clearly identified, corresponding to quasinormal modes with finite Q n . We 
will identify these peaks with quasiparticles. The different criteria existing in the 
literature to define a quasiparticle generally relate the imaginary part of the quasi- 
normal modes (r n , responsible of the width of the quasiparticle peaks) and energy 
(Q n related to the positions at which the peaks are centered), giving a condition of 
the form £^ < 1. Taking the T — > limit, these peaks can be seen to coincide 
with the supersymmetric mesonic spectrum (see for example [11] for the study of the 
transverse mode in the D3/D7 system). 

Ml = 2ir 2 M 2 n(n + 1), n>l, (4.4) 

with M the mass scale of the system. 

The quasinormal mode point of view turns out to be useful for understanding 
the qualitative behaviour of the spectral function in the quasiparticle regime under 
variation of the different parameters. We will give here some heuristic reasoning 
about these variations and compare it with numerical results obtained following the 
procedure described above. We will also link the behaviour of the quasinormal modes 
with the geometry of the D3/D7 system by means of the induced horizon on the 
probe branes. 

From figure 1 we can guess in what region of the parameter space the quasipar- 
ticle interpretation is appropriate. Large narrow peaks in the spectral function are 
associated with embeddings resembling Minkowski-like ones everywhere but in the 
region close to ip — 1> where a narrow throat, consisting of a bundle of fundamental 
Fl-strings pulling the brane into the horizon, is formed [10]. This narrow throat 
implies that the induced horizon on the probe brane has a very small area, corre- 
sponding to the flat region on the right of figure 1. From equation (3.15), we see that 
the quasiparticle regime is that of ipo « 1. In physical parameters this means that 
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the quark mass over temperature has to be high (and higher the larger the baryon 
density is). 

In the quasiparticle regime we expect equation (4.4) to roughly describe where 
the centers of the peaks of the spectral function should be. In this case one can 
identify the mass scale M with the mass of the constituent quarks by M = mT. 
Defining the mass of the (now melted) mesons using the dispersion relation of the 
quasinormal modes M 2 = u 2 — q 2 oc T 2 we observe that, when M increases, u: n has 
to increase correspondingly. 

Im[w] 

\ , , 1 - 1 - - ' — — — 1 — ' — Re[u>] 




Figure 4: Position of the quasinormal modes with positive real part as we vary m for 
fixed d = 0.01 and q = 0.01. The massless quark limit corresponds to the lower points on 
the graphs and we evolve up to m = 2.01. We see that when the quark mass is increased the 
pole gets closer to the real axis hardly changing the value o/O n . From a given value of the 
quark mass it changes completely the behaviour, approaching the real axis asymptotically in 
£l n . The very large frequency limit can be read as aT — > limit, therefore the poles should 
sit exactly on the real axis. The poles for different values of the parameters evolve in the 
same qualitative way. 



One expects that in the quasiparticle regime this qualitative behavior still re- 
mains valid, possibly with a slightly modified rate of change. This would mean that 
in the spectral function the peaks are shifted to larger values of the frequency, so the 
energy of the quasinormal modes, fl n , grows with increasing values of the parame- 
ter m = M/T. This is what we find numerically, as shown in figure 4 for a single 
pole. Notice also that an increasing value of M/T corresponds to a closer agreement 
with the quasiparticle condition ^f- <C 1 (see figure 1). This also supports the 
description given above in terms of the induced horizon, where at fixed T increasing 
the mass of the quarks M meant a smaller induced horizon, that is, the embedding 
of the probe branes resembles closely that of a meson in the non-deconfined phase 
when T/M — > 0. It should be noted that the low momentum modes which have 
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support over a large region of the D7-brane will see little effect from the narrow 
throat. Large spacetime momentum modes are concentrated around ip = I and will 
therefore notice the effects of the horizon even for small values of T/M. 

Another feature present in figure 4 is that, once we leave the quasiparticle regime 
of the theory, the position of the quasinormal modes evolve with the quark mass in a 
completely contrary way. The energy Q n associated to the quasinormal modes varies 
slightly, whereas the damping factor — T n is increased considerably. 

Im[uj] 




Re[u>] 



Figure 5: Position of the first eight quasinormal modes with positive real part as we 
vary d for fixed m = 2 and q = 3. The red points mark the values (beginning at the top) 
d = 0.01, 0.012, 0.026, 0.063, 0.135, 0.254, 0.432, 0.680, 1.01, 1.43, 1.96, 3.38, 5.37, 8.01, 
11.4, 15.6. Between any two consecutive red points there are ten data points. When d is 
increased all the quasinormal modes begin to orbit a certain point, but this happens beyond 
the quasiparticle regime, we have not investigated whether this is a numerical issue. 

The next parameter under consideration is the baryon density d. When the probe 
branes are charged the embeddings can be roughly described as being Minkowski- 
like with a throat entering the black hole. This way of seeing the embedding is more 
accurate the smaller d is (being associated to a narrower throat). Having only black 
hole embeddings, the mesons melt and we have a finite width for the peaks in the 
spectral function, corresponding to a finite value of T n . The peaks are broader the 
larger the induced horizon is. That is, when d is larger and the approximation of the 
embedding to a Minkowski-like embedding is worse. 

The conclusion is that the effect of increasing the baryon density on the quasinor- 
mal modes is to increase the value of | F n | , driving the system out of the quasiparticle 
regime. This is what we observe in figure 5. From the point of view of the induced 



- 28 - 



horizon area (figure 1 ) it is clear that an increasing baryon density d for a fixed value 
of m will broaden the peaks of the spectral function. 

Im[iu] 

, , j _ __, , , 1 , , , 1 , , , 1 — Re[ui] 




Figure 6: Position of some quasinormal modes with positive real part as we vary q for 
fixed m = 2 and d = 0.01. The red points correspond to q = 0.01 — 10.01 in regular steps. 
Between any two red points there are at least 10 data points. 

Now we turn our attention to the behaviour with q. The squared meson mass 
definition in equation (4.4) is M 2 = uj 2 n — q 2 . Therefore we see that if we want to 
keep M fixed as q increases, again the value of the frequency has to grow, meaning 
that Q n approximately grows with q. This is what we see in figure 6. 

As pointed out in [12], there is a maximum value of q at which the quasiparticle 
condition ceases to hold. In the cited paper the authors identify this critical value 
(\crit by studying the Schrodinger potential in the transverse vector channel. As an 
increasing q enhances the value of the energy of the quasinormal modes Q n following 
the dispersion relation, it is expected that before reaching q crit the value of the widths 
T n increases faster. In figure 7 we plot the continuation of figure 6 for higher values of 
the momentum. There we see how the modes enter a region where T n oc Q n , driving 
the system out of the quasiparticle regime (by diluting its effect on the spectral 
function). The value of the momentum at which this happens is different for each 
quasinormal mode. 

4.5 Dispersion relations 

Following the position of the quasinormal modes at large value of the momentum 
a study of the dispersion relations can be performed, giving information about the 
limiting velocity, v n , of the unstable quasiparticles associated to the modes. This is 
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Figure 7: Position of the quasinormal modes for m = 2 and d = 0.01. The red points 
correspond to q = 10.01 — 20.01 in regular steps. Between any two red points there are at 
least 10 data points. 

done by fitting the real part of these modes to a mass hyperbola 



An example of this can be found in figure 8. There we follow the first four quasinormal 
modes for m = 4 and d = 0.5. These parameters give sharp peaks in the spectral 
function, as can be guessed from figure 1. The numerical data can be fitted to a 
mass hyperbola like the one given above, which turns out to be a good fit. The 
dispersion relations in figure 8 correspond to the data shown. When q — > we can 
read the value of the mass of the quasiparticle produced. This is a decoupling limit, 
and one should expect that the quasinormal modes associated to the longitudinal 
electric field should give the same masses as the ones obtained from the transverse 
electric field, because of the recovery of rotational symmetry. This is indeed the case 
for the modes with masses M n = 4.02, 7.03 in the present example. 

Notice that the results given seem to suggest superluminal velocities, but the 
difference with the speed of light is so small that this effect cannot be differentiated 
from numeric instabilities reliably. A complete study of these dispersion relations 
would imply a great improvement of the numerical code, but this is beyond the 
scope of this work. 

4.6 Summary and conclusions 

We have described how to extract dynamical information from a system of coupled 
fields, which corresponds holographically to the description of operator mixing. This 




M 2 n + vW. 



(4.5) 
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Figure 8: Dispersion relation (a) and damping factors (b) for the first quasinormal modes 
as a function of momentum for m = 4 and d = 0.5. The limiting velocities are compatible 
with the speed of light. 



generalizes the celebrated (and widely used) prescription of [2] in a way that agrees 
with the causal structure of the retarded Green's function and with the general spirit 
of the AdS/CFT correspondence. We find that the matrix- valued Green's function 
for a coupled system is intimately related to the symmetries of the bilinear action 
describing the fields and their mixings. Given a generic action, taking fluctuations 
and linearizing the system gives rise to N 2 U(l) global symmetries (with iV the 
number of coupled fluctuations), the Green's function being given by the associated 
Noether current. 

As an example of the procedure we calculated numerically the quasinormal modes 
of the D3/D7 system in the presence of finite baryon density, where a coupling be- 
tween the longitudinal vector channel and the scalar channel appears. We described 
from the quasinormal mode paradigm how the system goes from a regime where the 
fields are independent to a regime where one has to consider the coupling between 
the fields, and how this is reflected in the spectral function. Understanding this we 
were able to calculate different physical properties, such as the diffusion constant 
and dispersion relations. 

It is also possible to see how the large time behaviour of the system makes a 
crossover between a behaviour dominated by a hydrodynamic mode and a quasi- 
particle mode. This transition is determined by a specific momentum q co (m, d) in 
the theory. When q < g co (m, d) the large-time behaviour of the system is set by a 
diffusive pole, whereas when q > q co (m, d) the large-time behaviour is described by 
the quasiparticles associated to the quasinormal modes. 

We were able to give a heuristic description of the behaviour of the spectral 
functions (and generally of the Green's function's non-analytic part) in the quasi- 
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particle regime, linking it to the behaviour of the quasinormal modes in this regime 
of parameters. We also gave a geometric point of view, relating it to the induced 
horizon of the flavor branes in the setup. It would be interesting to try to find phe- 
nomenological expressions relating the behaviour of the quasinormal modes in the 
quasiparticle regime and the induced horizon on the _D7-branes, but we have not 
investigated this here. 

The study of dispersion relations and the behaviour of the system at high fre- 
quency and momentum has still to be clarified. This matter has escaped the efforts 
of several groups due to the numerical instabilities that generically appear when us- 
ing numerical programmes with these range of parameters, ultimately related to the 
rapid oscillations of the system. 

We conclude stating that the general procedure introduced in this work is a very 
useful tool to treat the AdS/CFT models that continue to become more intricate 
as we simulate systems ever closer to the real world. In particular with the current 
interest in AdS/CMT models where the study of systems with coupled degrees of 
freedom is common such a procedure may help to understand a wealth of interesting 
phenomena. 
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A. General comments on the causal matrix Green's functions 

Consider the matrix-valued spectral function 



p tJ (x) = ([O l (x),O J (0)]), 



(A.l) 



which exhibits the following properties 



p( x y = p(- x ) = -p(x) t . 



(A.2) 
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Correspondingly, the Fourier transform p(k) = J d A k e tkx p(x) also satisfies a set of 
identities 

p(fc)t = ~ p {k) = -~ p {-kf . (A.3) 

In particular this means that the diagonal components are real and antisymmetric 
under k — > —k. One may also be interested in the behaviour under uj — > —uj. For 
0(3) invariant theories the diagonal components will also be real and odd in u 

pu(u, |q|) = pu(uJ, |q|)* = -pu(-U, |q|) . (A.4) 

For the off-diagonal components however, only if one also imposes time reversal or 
parity symmetry can one prove that the off-diagonal entries must be either even or 
odd functions of the frequency. In the present case time reversal symmetry is broken 
by the presence of a finite baryon density The parity operator acts as POi(t, x)P _1 = 
<TiOi(t, — x) with Oi = ±1, hence 

P[pij(t, X)] = (TiCTjPijit, -x) . (A. 5) 

Parity invariance implies pij(t,x.) = <Ji<jjPij(t, — x), which for the Fourier transform 
implies that 

p i:j (uj, q) = -(JiCJjPiji-UJ, q)* . (A. 6) 

So the off-diagonal entries are either odd or even functions of to depending on the 
signs Oi. In the case where the fields transform in the same way under the parity 
operator this means that the real (imaginary) part of the off-diagonal components is 
an odd (even) function of the frequency. 

From the spectral function, as defined in (A.I) we can define two causal propa- 
gators, namely the retarded and advanced Green's functions 

G R (x) = -i&(t)p(x) , (A.7a) 
G A (x) = iO(-t)p(x) , (A.7b) 

where x = (t, x). Using (A.3), one can prove the following relation amongst the 
Fourier transforms of these 

G R (k) = G R (-ky = G A (k)K (A.8) 

From here, we see that the real (imaginary) part, ReG^ (ImG^), is even (odd) 
under k — > —k. We can compute the Fourier transform of the retarded Green's 
function, which is given by the convolution of the Fourier transform of the Heaviside 
step function Q(uj) with the Fourier transform of the spectral function p(k), 

G R (u,q) = -i e(co-p)p(p,q)^. (A.9) 
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Using the Fourier transform of the step function 



oj + ie 

and the Sokhatsky-Weierstrass theorem we get 

G R (u, q )=V ^-^-—--p(u,q). (A.10) 

From the hermiticity of p(fc) we see that we can regard (A. 10) as a split of G R (k) 
into its hermitian and antihermitian parts, and find that the spectral function can 
be computed from the antihermitian part of the Fourier transform of the retarded 
Green's function 

p(k) = i[G R {k) - G R (ky] = 2zG { R A \k) , (A.ll) 

where the (A) stands for antihermitian 19 . Plugging this back into (A. 10) and taking 
the hermitian part (H) on both sides we arrive at 

G { n\co) = -vr^^dco', (A. 12) 

which is nothing but the Kramers-Kronig relation for the matrix Green's function. 
It is complemented by the conjugate relation interchanging the hermitian and anti- 
hermitian parts. Under parity transformation the Green's function satisfies 

G*{u^)=a l a 3 G? 3 {-u^y. (A.13) 



B. Formula for the residue 

Consider the adjugate matrix adjfif] defined by 

H- 1 = det^^adjftf] . (B.l) 

Note that adj[if] is finite at det[if] = 0. The relevant part of the Green's function 
is now 

d „\ 1 



G = 2A- -H . a dj[H] 



dr J L J det[#] 



F3.2 



which makes manifest that the poles are given by detfiJ] = 0. 

Close to a quasinormal mode, u n = Q n + iT n , the determinant as a function of 
the frequency can be approximated as 

d 

det[M (u;)] = (u - u n )—det[H] , (B.3) 



1!) 



Using (A. 8) we can always work with retarded Green's functions Gr. 
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since by definition det[H(uj n )] = 0. Therefore the matrix of residues is given by 

(B.4) 



TZ n = -2A-[ 4zH] ■ adjftf]^ 



dz J 1 J |-det[if] 



In practice one also faces the question of how to compute numerically the holo- 
morphic derivative Jj. The Cauchy-Riemann equations allow us to express the holo- 
morphic derivative of det[if (u;)] numerically as 

d , ,, m det[H(uj n + 5)}-det[H(u; n )} 

- det[H] = , (B.5) 

for a conveniently small and real S and where we have not taken det[H(cj qnm )} = 
because it is numerically more stable. 



C. Boundary action 

Expanding the BI lagrangian (3.1) up to second order in fluctuations, and performing 
a Fourier transformation on them as in (2.2) we may cast the result in the form 
given in (2.5) for the gauge invariant fluctuations (E^(u),^k(u)) and E%(u) which 
decouple. From here we can read off the coefficient matrices Au, Bjj and Cu that 
can be seen in (2.5). Notice however that with this lagrangian, the equations of 
motion lead to an asymptotic expansion for ^(u) akin to the one given in (3.11), 

or *$?k( u ) U —>° clt/u + In the spirit of the discussion at the end of section 2.1 we 

must rescale the fluctuations by a matrix $ 7 = D 1 j§ J , or 



E£(u)\ = /lfl \ (E£(u) 



(C.l) 



and compute the new coefficient matrices Au, Bjj and Cu as given in (2.8a), (2.8b) 
and (2.8c). The coefficients, Fjj, are functions of the embedding solution, d, u, q 
and u which are best given in terms of the following functions 



h(u) = ^ ^ ^ . (C.2b) 



g(u) = ^){uf + d 2 u 3 , (C.2a) 

#H 

We find, in terms of the usual dimensionless ratios ro = u>/2irT and q = qjl-nT 
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LL ^-q 2 /^ 6 (27rT)2 u * - q 2 + 

^iqduWphW 10 (itTL 2 \ 1 ^ 9 . , n a g /ttTL 2 \ 



with 



- q 2 /^ e 

, . r fiP M h 3 / 2 (tt> 2 iP 2 g - q 2 (fi> 8 - 4d 2 /VV)) /ttTL 2 ^ 2 

Axpiif — TV ~ — — - 

g(^g - q 2 /V 6 ) V 2 ™' 

~ A,(^) 2 + ..., (C.3, 

-Bll = B L q, = , 

-iqd^Ifi> 8 h(2ufi> 2 g^' + ^(3^ 8 + Au 2 f{d 2 u 3 + 4^ 6 )^' 2 ) [kTL 2 \ 1 



g(-q 2 ft/j 6 + rc 2 #) V 2 ™' / 27rT 

/^ 12 h 3 / 2 ( <kTL 2 \ 2 
g(-q 2 f& + K> 2 g)u \ 2na' J 

ra 2 (# 4 + 2u{2u 3 d 2 - i> 6 )i> 2 ipip' + 8u 3 f(d 2 u 3 - 2^ 6 )^' 3 ) 

+ q 2 /(-^ 8 (^ 2 - 2mM)') + 4u 5 d 2 f^' 2 + 8u 3 f(d 2 u 3 + 2^ 6 )^' 3 )) , 



AT AT rp2 

M = N f T D7 Vol(S 3 )(2Tra') 2 (irTL 2 ) 2 = f c . (C.4) 



D. Equations of motion for the fluctuations 

In this appendix we reproduce for completeness the equations of motion for the 
fluctuations and analyze with care the limit q — > 0. As was shown in [13], the set of 
fluctuating fields Aq,A\ and \P can be shown to satisfy a set of coupled differential 
equations for the gauge invariant combination (longitudinal electric field) El given 
in (3.18) 

El + %E' L + £E L + + ' + m> = , (D. la) 

V" + SV + ®V + 4E' L = 0. (D.lb) 



with 



% = % ll » = «i ; £ = Ci ; 56 = *)!; C = d ; (D.2a) 



^C^^C^^C^' (D ' 2b) 
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where 



*2 

Ci 
€ 2 

£> 2 



log' 



u 2 + q 2 f 



log' [v^rr^] 



a ,2 7 oo 



a; 2 7 00 + 5 2 7" 



+ 



A'-S 



q 2 Y 



A cj 2 7 00 + g 2 7" 



^2 



^2 7 00 + g 2 7 n 



r 



97 



Or 



A 



g(l-^'A) 



97 



^'7 0r 

Or 



^OOyr 



(S + A) 



„0r 



7 A log' ' 7 



k^^UU _|_ q^ryll ryVT - I ^ 

9(1 " t A) log' [V=T7"-^(1 - ^'A)] 



oo 



w 2 7 00 



57 



Or 



^2^00 _|_ g2 7 « ^OO^rr 



(log' (v/=77 il 7 rr ) - 



«i 

«2 



97 



Or 



^OO^rr 



97 



Or 



^OO^rr 



S'- A 

00 



cu 2 7 00 + gV 



qu 2 7 0r _ , / 7 

;S log 



2 — A 

ryrr /j 2 -y00 _|_ ^2^11 



^rr / CJ 2 7 00 + g 2 7** 7 rr \7°° 

,S log'^V) 



q(l - if/ A) 



H(z), 



where 



1 



E = - ( 7 U V 2 GW - 3 7 nQ G'nn,v-) 5 A = 7 "VGW , 



with G^y components of the original 10-dimensional bulk metric and H(u) is a rather 
lengthy expression we give here for completeness 



H{u) 



(cu 2 7 00 (l - 



J/^GW^ i „2„,11 



-77 M "C^(1 - V'A) 

•n '''A)) 



7 UU (1 - V>'A) 

(l (7° n GW) 2 + Wt^^G^Gmw + 3 7 nn G nn ,^) 

27" M G^(1 - V'A) 
[ip' 2 G^ t ^ — \^ uu (G^^) 2 ip' 4 ^ 
2G^(1 - iSA) • 
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Written in this form the equations (D.la) and (D.lb) decouple smoothly in the 
limit q — > 0. We have used the following definitions for the background matrix 20 
lab = g a b + 27ra'F ab 



2ira'A' (u) , 7' 



00 



uu 



7oo 




iL + 7oo7 
1 



uu 



7 



70n + 7oo7 
-7o« 



uu 



7 



On 



7o« + 7oo7 



uu 



Note also that 7o« = — 7«o and 7°" = — 7"°. We also denote ^—7 = v /— det 7 a fe 



taking into account only the radial part. 
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